Width Distributions and the Upper Critical Dimension of KPZ Interfaces 
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Simulations of restricted solid-on-solid growth models are used to build the width-distributions 
of d — 2 — 5 dimensional KPZ interfaces. We find that the universal scaling function associated with 
the steady-state width-distribution changes smoothly as d is increased, thus strongly suggesting that 
d — 4 is not an upper critical dimension for the KPZ equation. The dimensional trends observed in 
the scaling functions indicate that the upper critical dimension is at infinity. 
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The KPZ equation has been introduced to model 
growth in terms of a moving interface. The equation is 
written for the height h(r, t) of the interface above a d- 
dimensional substrate 

d t h = i>V 2 h + \(Vh) 2 +r], (1) 

where v and A are parameters, while rj(r, t) is a Gaussian 
white noise. Eq. (|l|) can also give account of a number 
of other interesting phenomena (Burgers turbulence, di- 
rected polymers in random media, etc.) and, accordingly, 
a lot of efforts has been spent on finding and understand- 
ing the scaling properties of its solutions ■ These in- 
tensive studies notwidthstanding, a number of unsolved 
issues remain, the question of upper critical dimension 
{d u ) being the most controversial one. 

The importance attached to d u stems from the hope 
that, in analogy with equilibrium critical phenomena, a 
better understanding can be achieved through system- 
atic expansions in terms of d u — d. The search for d u has 
been on for about a decade f5j-|l6|l and the results range 
from d u « 2.8 to d u — oo. Analytical estimates originate 
mainly from mode coupling theories which yield exact 
results for d — 1 fll7|| . Extending this approach to higher 
dimensions |7|-|ll| one obtains values of d u which, after 
refining the self-consistency schemes, appear to settle to 
d u = 4. The result d u = 4 also emerges from various 
phenomenological field-theoretic schemes p,[L2| and some 
nontrivial consequences of the phenomenological argu- 
ments appear to be in agreement with simulations Eq] . 

In contrast to the analytical approaches, numerical so- 
lution of the KPZ equation GJ|, simulations of systems 
belonging to the the KPZ universality class |l4] , |l5| ] , and 
the results of real-space renormalization group calcula- 
tions pi provide no evidence for a finite d u . Further- 
more, the only numerical study Eq ] of the mode-coupling 
equations gives no indication for the existence of a finite 
d u either. 

There are, of course, problems with both the analytic 
approaches and the numerical works. Assumptions about 



the scaling structure of the solution underlie the field 
theoretic approaches, and uncontrolled approximations 
are made when writing down the governing equation in 
mode-coupling theories. Additional uncertainties come 
from the use of various selfconsitency schemes in solving 
the mode-coupling equations. Simulations and numerical 
works have their own share of difficulties. The systems in 
higher dimensions cannot be large; the extraction of ex- 
ponents using fitting procedures which involve unknown 
correction-to-scaling exponents makes the error estimates 
suspect, and there may be difficulties with the numerical 
solution of the mode-coupling equations as well O ], 

In view of the above controversy, it is highly desirable 
to approach the d u problem in a way unbiased by ap- 
proximations and fitting procedures. Such an approach is 
described below where we study the steady-state width- 
distributions of d = 1 to d = 5 KPZ interfaces. 

The width-distributions have been introduced to pro- 
vide a more detailed characterization of surface growth 
processes |l9|-|22]|, and they have been used to establish 
universality classes of rather divers phenomena ]^3| p8[ . 
The quantity whose distribution is of interest here is the 
mean-square fluctuation of the interface defined by 

W2 = ±-Y<[h(r,t)-h} 2 , (2) 

r 

where Al is the area of the substrate of characteristic lin- 
ear dimension L, and h = '^ l fh{r,i)/Ai, is the average 
height of the surface. Sampling u>2 in the steady state, 
one can build the so called width distribution Pr,(w2)dw2, 
defined as the probability that W2 is in the interval 
[u>2, w 2 + dw 2 ]. If the quantities h(r, t) — h were uncor- 
rected at large distances the probability distribution of 
W2 would be approaching S(w2 — (w2)l) for L — » oo. On 
the contrary, the fact that the distribution is non trivial 
implies that these quantities are strongly correlated at 
large distance. 
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The usefulness of this distribution lies in the follow- 
ing observation supported by all the examples studied 
so far (including d — 1 and 2-dimensional KPZ surfaces 
]Td| ^l],|7)). Namely, in systems where the steady-state 
roughness diverges (w2)l — * oo in the L — > oo limit, 
Pl{w2) assumes a scaling form 



Pl(u> 2 ) 



1 



(w 2 ). 



W2 



(3) 



where §d(x) is an universal scaling function characteris- 
tic of the universality class of a given nonequilibrium dy- 
namics in dimension d. This universality is understand- 
able, it is a consequence of the facts that (i) a steady 
state can be considered as a critical state if the fluctua- 
tions diverge, and (ii) in critical systems, the distribution 
functions of macroscopic quantities (such as (1U2)) are 
characterized by scaling functions which are universal. 

The universality of &d(%) allows the investigation of 
the problem of d u , once it is noted that the scaling func- 
tions depend on dimensionality up to d = d u and they are 
expected to take on a fixed shape for d > d u - Thus if one 
finds that scaling functions vary smoothly in dimensions 
1 < d < d, one can conclude that d — 1 < d u . This is 
the line of argument we employ below for KPZ systems. 
We shall compare the scaling functions (||) for 1 < d < 5 
using the exact results for the d = 1 steady state [ p^| , 
previously obtained simulation data for 2 < d < 4 re- 
stricted solid-on-solid (RSOS) growth models and 
by generating new data for the d = 5 RSOS model. Our 
main finding is that the $d(x)-s change smoothly as d is 
varied thus suggesting that d u > 4. 

It is important to recognize that there are no fit- 
ting procedures in the above approach. The width- 
distributions are just histograms calculated from Monte- 
Carlo (MCI simulations. Both quantities W2 and (11)2) L 
entering (|3j) are measured and no scaling properties of 
{1V2) L are used or assumed. The only approximation 
is the finite size of the systems investigated. It should 
be noted, however, that our approach relies only on the 
shape of the scaling functions. Since the important size- 
dependences reside in the argument of these functions the 
functional forms converge at small sizes. A further and 
rather important observation that helps to reach our con- 
clusion is that, as we shall show below, the scaling func- 
tions converge to well distinguished forms for 1 < d < 5. 

Let us now present and discuss the evidence for our 
conclusion of d u > 4. The scaling functions (|^) for di- 
mensions d = 1 — 5 are displayed in Figs. 1 and 2. The 
d = 1 curve is an exact result The rest is obtained 
from simulations of the RSOS model that is believed to 
belong to the universality class of the KPZ equation p9| . 
The RSOS model and its simulations are described in 
Jl5| where hypercubes of volume L d with periodic bound- 
ary conditions were simulated and a multi-surface coding 
technique allowed to obtain excellent steady-state statis- 
tics for systems up d = 4. We take the results from this 



work to build $>d{%) for e? = 2 — 4 and, in this paper, we 
extend these simulations to find $>d( x ) for d = 5 as well. 
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FIG. 1. Scaling function for the width-distribution [eq. (|3|)] 
of KPZ systems in dimensions d = 1 — 5, with L giving the 
largest linear size of the system in which <&d(x) was measured. 

As one can see on Figs. 1 and 2, the scaling functions 
change smoothly as d increases. The 4>d(x)-s get nar- 
rower and more centered on x = 1, and there does not 
seem to be any break in this behavior at d = 4. The 
equality of the d=4 and 5 scaling functions appears to 
be excluded. Since our conclusion about d u > 4 rests on 
the above observations we must now discuss some details 
in order to make it more than a visual observation. 
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FIG. 2. The same as Fig.l but on a semilog scale in order 
to demonstrate the differences at small probabilities. Note 
that the range of x has been enlargened so that the large- a; 
asymptotics would be better seen. 

The basic problems that may arise in measuring steady 
state properties are the problems of statistics, relaxation, 
and finite-size. Since the multi-surface coding allowed 
the simulations of 32 or 64 systems in one run, we had 



2 



no problem gathering data with good statistics. The re- 
laxation time problems were taken care by having very 
long runs and being in the asymptotic plateau region of 
u>2 for at least an order of magnitude longer period than 
the time of reaching the plateau (for details see discussion 
and Figs. 3, 5, and 7 in 

The solution to the problem of finite size is less obvi- 
ous. In general, one can observe that the $d-s converge to 
their limiting shape when the number of sites (N = L d ) 
becomes about a factor 2 larger than the number of sites 
on the surface (N a ). Fig.|] demonstrates this observation 
for the d = 2 and 4 systems. Results for linear sizes of 
L = 7 (N w 2N a ) and L = 157 (N w A0N s ) are com- 
pared for d = 2 and one can see that they have the same 
$-s within the the statistical errors of the simulations 
fjof . Similar conclusion can be drawn by comparing the 
L = 13 (N w 2N S ) and L = 25 (N « 3.5/Y s ) systems in 
d = 4 (analogous results for d = 3 are not displayed on 
Fig.^ in order to keep clarity in presentation). 
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FIG. 3. Finite-size effects on the scaling function [eq.(|^)] 
in dimensions d = 2, 4 and 5. L denotes the linear size of the 
hypercubes investigated. 



We need the N « 2N S convergence rule because the 
largest d = 5 system we can investigate has L = 15, cor- 
responding to TV w 2N S . The results for <j> 5 displayed in 
Fig.|| indicate that the N ss 2N S rule applies to d = 5 as 
well. Indeed, systematic deviations between the L = 11 
and L = 15 curves can be detected only at small values 
of $5 in the region of x < 0.85. An important feature 
of the size-dependence of $5 that can be seen in Fig.|| 
is that the maximum of $5 slightly increases with size. 
This means that, near the maximum, $5 — $4 becomes 
larger with increasing L, thus excluding the possibility of 
the two functions $5 and $4 becoming equal. The dif- 
ferent functional forms for the scaling functions in d = 4 
and 5 then indicate that d = 4 is not the upper critical 
dimension for the KPZ systems. 

The concept of smooth changes across d — 4 can be put 



on a more quantitative basis by examining the dimen- 
sional trend in the spread of the scaling function around 
its average x = 1 

o\= [™ dx(x-l) 2 ^x) = ^±-l 7 (4) 
Jo {W2) 1 

that is related to relative mean-square fluctuations of W2 ■ 
Apart from the case of d = 1, the L dependence of a c i 
is very weak and plotting Ud{L) against l/L one can get 
accurate estimates of <7,i(oo). The values of 0-^(00) as a 
function of X/d are displayed on Fig.f§ As one can see, 
apart from the d = 1 result, the straight line ad ss 0.71/d 
gives an excellent description of the dimensional depen- 
dence of ad- 




1/d 

FIG. 4. Dimensional dependence of the relative fluctua- 
tions of u>2- The extrapolated values Od{L — > 00) — 0^(00) 
[see eq.(Q)] are plotted against 1/d with the solid line given 
by a d = 0.71/d. 

The result ad ~ 0.71/d indicates that ad —> for 
d — > 00 i.e. the scaling function converges to a delta- 
function, 8(x), at d = 00. Remarkably, the convergence 
&d(x) — > S(x) also takes place in a related surface- growth 
model, in the Edwards- Wilkinson model if d — > 2 [ po[ . 
Since d = 2 happens to be the upper critical dimension 
of this model (the interface becomes flat for d > d u = 2) 
one may speculate that the results displayed on Fig.H 
actually give support to the suggestion that failure of 
numerical attempts at locating a finite d u means that 
d u = 00 for KPZ systems |pl| . 

In summary, we believe that the main results of this 
paper (Fig|l] and |^) provides strong evidence that d u > 5 
and, furthermore, Fig.^ suggests that d u = 00. As a final 
remark, let us note that the results displayed on Fig.^ and 
^ can be appreciated from another point of view. Namely, 
we have constructed the scaling functions of width distri- 
butions for the KPZ universality class. Thus we have ex- 
panded the picture gallery of scaling functions that may 
be used for identifying the universality classes of nonequi- 
librium steady states. 
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